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Abstract 

Background: Brucellosis caused by Brucella abortus is one of the most important zoonoses in the world. 
Multiple-locus variable-number tandem repeat analysis (MLVA16) has been shown be a useful tool to epidemiological 
traceback studies in B. abortus infection. Thus, the present study aimed (/') to evaluate the genetic diversity of B. abortus 
isolates from a brucellosis outbreak, and (/'/) to investigate the in vivo stability of the MLVA16 markers. 

Results: Three-hundred and seventy-five clinical samples, including 275 vaginal swabs and 1 00 milk samples, were 
cultured from a brucellosis outbreak in a cattle herd, which adopted RB51 vaccination and test-and-slaughter 
policies. Thirty-seven B. abortus isolates were obtained, eight from milk and twenty-nine from post-partum/abortion 
vaginal swabs, which were submitted to biotyping and genotyping by MLVA16. Twelve B. abortus isolates obtained 
from vaginal swabs were identified as RB51. Twenty four isolates, seven obtained from milk samples and seventeen 
from vaginal swabs, were identified as B. abortus biovar 3, while one isolate from vaginal swabs was identified as 
B. abortus biovar 1. Three distinct genotypes were observed during the brucellosis outbreak: RB observed in all isolates 
identified as RB51; W observed in all B. abortus biovar 3 isolates; and Z observed in the single B. abortus biovar 1 isolate. 
Epidemiological and molecular data show that the B. abortus biovar 1 genotype Z strain is not related to the B. abortus 
biovar 3 genotype W isolates, and represents a new introduction B. abortus during the outbreak. 

Conclusions: The results of the present study on typing of multiple clinical B. abortus isolates from the same outbreak 
over a sixteen month period indicate the in vivo stability of MLVA16 markers, a low genetic diversity among B. abortus 
isolates and the usefulness of MLVA16 for epidemiological studies of bovine brucellosis. 
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Background of sanitary barriers to international trade of animals and 

Brucellosis is a bacterial disease caused by microorgan- their products [1]. 

isms of the genus Brucella that affects different live- Because direct and indirect losses triggered by B. abortus 
stock and wild animal species, besides man [1]. Cattle infection, the control and eradication of bovine brucellosis 
are the preferred host of Brucella abortus and the eco- is an important goal of several countries where the disease 
nomic importance attributed to bovine brucellosis is is endemic, including Brazil that since 2001 has a national 
based on losses caused by abortions, stillbirths, weight program on the control and eradication of brucellosis and 
loss, decreased milk production and the establishment tuberculosis (Programa National de Controle e Erradicacao 

da Brucelose e Tuberculose - PNCEBT) [2,3] . 
Typing of Brucella spp. by biovar determination and 
* Correspondence: alageiavet.ufmg.br genotyping are important tools in a brucellosis control 

'Laboratorio de Bacteriologia Aplicada, Departamento de Medicina and eradication programs [4]. Indeed, molecular typing 

; 11 ! ; ^ 1 1 ' " 1 r methods are commonly used to investigate epidemio- 

logical relationships among isolates and sources of 
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infection [5]. In this sense, the multiple locus variable 
number of tandem repeats (VNTR) analysis (MLVA) has 
proved to be an important tool in molecular epidemiology 
studies of brucellosis [6-11], in the characterization of new 
Brucella species [12] and in the evaluation of sources of 
human infection [13]. 

In MLVA context, the assessment of the genetic stabil- 
ity is one of the vital elements to guarantee the success- 
ful use of this typing method, since it could confirm its 
ability to differentiate between two independent strains, 
while depicting relationships among strains from the 
same source [14]. The genetic stability can be deduced 
by typing results of multiple isolates originating from the 
same strain or strains after numerous passages [8,9,15,16]. 
In in vitro evaluations, B. abortus vaccine and field strains 
have been shown to be stable by MLVA16 [9,15-17]. 
Regarding to in vivo stability of the MLVA16 markers, 
data based on B. abortus experimental infection showed 
only few changes in the hyper- variable loci [9]. Whereas, 
MLVA16 in vivo stability have not been evaluated under 
natural conditions. 

Therefore, the aims of this study were (/) to evaluate 
the genetic diversity of B. abortus isolates obtained from 
a single cattle brucellosis outbreak and (u) to investigate 
the in vivo stability of the MLVA16 markers. 

Results 

Serologic testing 

The incidence of brucellosis in the herd between May 2009 
and January 2011 is shown in Figure 1. No positive animal 
was observed on three monthly tests after October 2011. 



Isolation and identification 

Thirty-seven isolates of B. abortus were obtained from 
340 cows (375 samples) during the 16 months period of 
outbreak: eight from milk samples (8% - 8/100) and 
twenty-nine from post-partum vaginal swabs (10.5% - 
29/275) (Table 1). Of the B. abortus isolates obtained from 
vaginal swabs samples eight were associated to abortion/ 
stillborn occurrence. Among animals that had vaginal 
swabs and milk samples collected (35 animals) three B. 
abortus isolates were obtained, all from vaginal swabs. 

Based on phenotypic characteristics, twenty four strains 
were classified as B. abortus biovar 3 (seven from milk and 
seventeen from vaginal swabs) and thirteen as B. abortus 
biovar 1 (one from milk and twelve from vaginal swabs) 
(Table 1). All strains classified as B. abortus biovar 3 pre- 
sented smooth colonial morphology. Among the thirteen 
B. abortus biovar 1 strains, twelve exhibited rough colonial 
morphology and growth on tryptose agar plates supple- 
mented with rifampicin (200 ug/mL) (all from vaginal 
swabs), whereas one strain isolated from a milk sample 
was classified as smooth Brucella spp. 

All strains were confirmed as belonging to the genus 
Brucella by genus-specific PCR [18], and as B. abortus by 
AMOS-enhanced PCR and Bruce-ladder multiplex PCR 
[19,20]. Moreover, the twelve isolates phenotypically clas- 
sified as B. abortus biovar 1 rough strains were confirmed 
as RB51 vaccine strain by AMOS-enhanced PCR [19]. 

MLVA16 

The set of genetic markers which comprises the MLVA16 
was observed to be stable in vivo, even when evaluated 




Figure 1 Monthly incidence of serologically positive cattle in a bovine brucellosis outbreak from May 2009 to January 201 1, Minas 
Gerais, Brazil. All cattle older than 24 months were monthly tested for brucellosis [Rose Bengal Plate Agglutination Test (RBPAT) as a screening 
test and the Standard Tube Agglutination Test (STAT) and 2-Mercaptoethanol Test (2ME) as a confirmatory test]. Positive animals were culled 
until three consecutive negative herd results were obtained. The herd was vaccinated with RB51 on May 2009 and on December 2009 (RB51 
vaccinations are highlighted with arrows). 
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Table 1 Clinical signs, biotyping, 
brucellosis outbreak 



molecular identification and genotyping of B. abortus isolates from a bovine 



Strain 


Clinical sample 


Collection date 


Biovar 


Genoty 


1414 


vagina 


swab 


09/01/2009 


3 


W 


2757 


vagina 


swab 


09/01/2009 


1 


RB 


2985 


vagina 


swab 


09/01/2009 


3 


W 


3144 


vagina 


swab 


10/01/2009 


1 


RB 


2575 


vagina 


swab 


09/02/2009 


3 


W 


2558 


vagina 


swab 


07/03/2009 


1 


RB 


2608 


vagina 


swab 


08/03/2009 


1 


RB 


1379 


vagina 


swab 


07/04/2009 


1 


RB 


3136 


vagina 


swab 


1 0/05/2009 


1 


RB 


2634 


vagina 


swab 


1 0/07/2009 


3 


W 


1365 


milk 


07/27/2009 


3 


W 


2476 


milk 


07/27/2009 


3 


w 


2598 


milk 


07/27/2009 


3 


w 


2683 


milk 


07/27/2009 


3 


w 


2785 


milk 


07/27/2009 


3 


w 


2965 


milk 


07/27/2009 


3 


w 


2980 


milk 


07/27/2009 


3 


w 


3159 


vagina 


swab 


07/27/2009 


1 


RB 


2922 


vagina 


swab 


08/1 8/2009 


1 


RB 


3169 


vagina 


swab 


08/21/2009 


1 


RB 


2545 


vagina 


swab 


09/1 9/2009 


3 


W 


2595 


vagina 


swab 


10/17/2009 


3 


w 


1198 


vagina 


swab 


1 0/28/2009 


3 


w 


1219 


vagina 


swab 


1 0/28/2009 


3 


w 


2251 


vagina 


swab 


07/11/2009 


3 


w 


2496 


vagina 


swab 


07/11/2009 


3 


w 


2700 


vagina 


swab 


09/11/2009 


3 


w 


2588 


vagina 


swab 


01/28/2010 


3 


w 


3238 


vagina 


swab 


05/03/2010 


3 


w 


2726 


vagina 


swab 


03/28/2010 


1 


RB 


3006 


vagina 


swab 


05/04/2010 


3 


w 


2998 


vagina 


swab 


04/15/2010 


3 


w 


3022 


vagina 


swab 


04/27/2010 


3 


w 


3016 


vagina 


swab 


04/30/2010 


3 


w 


2774 


vagina 


swab 


02/06/2010 


1 a 


z 


2543 


vagina 


swab 


03/11/2010 


I 


RB 


2850 


vagina 


swab 


03/11/2010 


1 


RB 



Clinical signal 

Abortion 



Abortion 
Stillbirth 

Stillbirth 



Stillbirth 

Abortion 
Abortion 

Abortion 
Stillbirth 



a 6. abortus biovar 1 field strain. 
b MLVA16 genotype. 

from different clinical specimens over a 16 months out- 
break. Analysis of the MLVA16 loci revealed three distinct 
genotypes among the 37 B. abortus strains isolated during 
the brucellosis outbreak, herein these genotypic pat- 
terns were labeled RB, W and Z, for new genotypes the 



identification adopted was sequential in the same way 
used by Minharro et al. [11] (Table 1). The genotype RB 
was observed in all twelve strains identified as B. abortus 
RB51, which were all obtained from vaginal swab samples. 
The genotype W was observed in all twenty four isolates 
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classified as B. abortus biovar 3 (seven from milk and 
seventeen from vaginal swabs), while genotypic profile Z 
was observed in only one strain, which was isolated from 
milk and identified as a smooth B. abortus biovar 1 field 
strain. This biovar 1 strain was isolated eleven months 
after confirmation of brucellosis in the herd on February/ 
2010, being one of the last strains isolated, since clinical 
samples were collected until May/2010 (Table 1). 

Patterns obtained in the sixteen VNTR loci are sum- 
marized in Table 2. Comparison of the results observed 
in the eight conserved loci included in panel 1 by Le 
Fleche et al. [6] with those available in the MLVA bank 
2014 (http://mlva.u-psud.fr/brucella/) showed genotypes 
27, 40 and 28 for the strains classified as MLVA 16 
genotypes RB, W and Z, respectively. Genotype RB was 
identical to the previously describe [6,16] for RB51 
vaccine strain. The MLVA16 pattern for genotypes 
W and Z did not find correspondence with those depos- 
ited on MLVA bank 2014. In comparison with genotype 
W, the genotype Z exhibited different alleles in the 
loci Bruce04, Bruce06, Brucell, Brucel6, Brucel8 and 
Bruce30 (Table 2). 

Discussion 

Results from the present study on 37 B. abortus isolates 
from milk and vaginal swabs taken from 340 cows (375 
samples) during a 16 months period showed the in vivo 
genetic stability of the MLVA16 markers. Molecular 



Table 2 Allelic types of MLVA16 loci observed for 
B. abortus strains isolated from a bovine brucellosis 
outbreak 



Locus 


Tandem repeat copy numbers 


Number 
of alleles 


Panel 


Genotype W 


Genotype Z 


Genotype RB 


Bruce06 


3 


4 


4 


2 




Bruce08 


5 


5 


5 


1 




Brucell 


3 


A 


A 


2 




Bruce12 


12 


12 


12 


1 




Bruce42 


2 


2 


2 


1 




Bruce43 


2 


2 


3 


2 




Bruce45 


3 


3 


3 


1 




Bruce55 


3 


3 


3 


I 




Bruce18 


/ 


6 


6 


2 


2A 


Bruce19 


42 


42 


42 


I 


2A 


Bruce21 


8 


8 


8 


I 


2A 


Bruce04 


5 


3 


3 


2 


2B 


Bruce07 


4 


4 


/ 


2 


2B 


Bruce09 


3 


3 


3 


1 


2B 


Bruce 16 


3 


5 


3 


2 


2B 


Bruce30 


3 


6 


5 


3 


2B 



typing methods are commonly used to investigate epi- 
demiological relationships among isolates and sources 
of infection [5], However, before being used for those 
purposes, PCR methods for molecular typing require 
careful in-house validation of typeability, reproducibility, 
repeatability, stability, discriminatory power and epide- 
miologic concordance [5,21]. The findings of this study, 
associated with previous data on the high discriminatory 
power and epidemiologic concordance of MLVA16, 
besides its good typeability and in vitro stability 
[6,7,9-11,13,15-17], corroborate the use of MLVA16 as 
suitable typing method for refining the understanding of 
the epidemiology of bovine brucellosis. 

Our results showed a low genetic diversity and the 
existence of three different B. abortus strains within a 
focus of bovine brucellosis, the RB51 vaccine strain 
(genotype RB), recovered from animals that were twice 
vaccinated in the period, and two field strains (geno- 
types W and Z). 

In the present study no RB51 was isolated from the 
milk of any animal during the post-partum period 
(30-60) but only from vaginal swabs, although some 
cows have been vaccinated in the last third of preg- 
nancy and a highly sensitive diagnostic strategy have 
been employed [22]. These results were corroborated 
by previous studies that also reported no recovery of 
B. abortus RB51 from any of the milk sample tested by 
conventional bacteriological methods of cows vaccinated 
during pregnancy or 30-60 days after delivery [23,24]. 
Regarding the presence of viable RB51 in postpartum 
vaginal secretion, it is important to consider that all ani- 
mals from which RB51 were recovered after delivery 
were over 5 months of gestation by the time of revaccin- 
ation with this strain and that no abortion was related 
to RB51 isolation (Table 1). Thus, it is likely that the 
RB51 booster had contributed to the recovery of the 
vaccine strain in postpartum period, since it has been 
observed that during transition period a depression in 
cell-mediated response occurs, which leads to a de- 
crease of the resistance to disease or increase of the re- 
sidual virulence of vaccines [25]. 

Other important finding of this study was the demon- 
stration of viable B. abortus in milk from infected ani- 
mals associated with the outbreak, corroborating the 
public health risk of the consumption of raw milk and 
unpasteurized dairy products [26]. The colonization of 
the mammary gland and associated lymph nodes by B. 
abortus with excretion of microorganisms in milk was 
already demonstrated [27]. In addition, it is important to 
emphasize that the methodologies used for bacterial cul- 
ture and molecular identification were able to differenti- 
ate RB51 vaccine strain from B. abortus field strains. 

Interestingly, field strains isolated from this brucel- 
losis outbreak showed distinct genotypes and biotypes, 
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suggesting that the outbreak had two different sources 
of infection. The majority of B. abortus isolates from 
the outbreak (24/25) were classified as B. abortus biovar 
3 genotype W, which remained the single cause of bru- 
cellosis in the herd for eleven months (from March 
2009 to February 2010). This widely demonstrates that 
the introduction of this strain in the herd was respon- 
sible for the occurrence of the outbreak and conse- 
quently for the high abortion rate observed at the end 
of 2008. Furthermore, in February 2010, the introduc- 
tion into the herd of another field strain, B. abortus 
biovar 1 genotype Z, was observed. This genotype Z was 
confirmed by a twice repetition of MLVA16 genotyping 
assay that showed the same results. There are two possible 
explanations for the presence of this infected animal, 
negative to the conventional serological tests adopted in 
the herd: first, this heifer was congenitally and persistently 
infected without seroconversion [28], and second, the 
negative serological result from this infected heifer was 
a false negative result inherent to any diagnostic test. 
However, the diagnostic strategy employed was very 
sensitive [29], which tends to minimized false negative 
results. Since this second field strain was represented by 
just one isolate, it is very likely that the control policy 
adopted in the herd, which included test-and-removal 
procedures, has prevented the dissemination of this 
strain to other animals and consequently a new increase 
in the brucellosis prevalence levels into the herd. 

Besides the differences in the classification of biovar, 
the differences observed between the genotypes of the 
field strains were very large and not limited to differences 
in the more variable loci (Bruce04, Brucel6, Brucel8 and 
Bruce30), but were also observed in conserved loci, such 
as Bruce06 and Brucell (Table 2). Therefore, the identi- 
fication of the new source of infection in the herd was 
only possible due the use of MLVA16, otherwise the 
genotype Z strain would have been considered just an- 
other B. abortus isolate within the outbreak. 

Epidemiological data of the herd also confirmed that 
the B abortus biovar 1 genotype Z strain was a newly in- 
troduced strain in the farm, since the young heifer from 
which this strain was isolated was introduced in the herd 
shortly at the end of the outbreak (November/2009). In 
fact, it has been widely demonstrated that the purchase 
of infected animals is the main risk factor for the intro- 
duction of brucellosis in free herds [30]. Unfortunately, 
it was not possible to obtain more epidemiological data 
about the heifer from which the second B. abortus field 
strain was isolated, such as the brucellosis status of 
the herd of origin, which would have allowed a better 
understanding about the epidemiology of the genotype 
Z strain. 

The MLVA16 panel 1 profiles of both B. abortus iso- 
lates, genotypes 28 and 40 by the MLVAbank, have 



already been previously observed in Minas Gerais State 
[11]. Moreover, the comparison between MLVA16 geno- 
types of B. abortus biovar 3 field strain (W) and those 
previously described by Minharro et al. [11] revealed 
differences restricted to hypervariable loci [panels 2A 
(Brucel9) and 2B (Bruce04)], suggesting a possible epi- 
demiological link between these strains, since all were 
isolated in the state of Minas Gerais and were also iden- 
tified as biovar 3. For the single field strain classified as 
genotype Z, the comparison with MLVA16 patterns 
previously described by Minharro et al. [11] showed a 
genetic distance of one locus (Brucel9) from a B. abortus 
biovar 2 also isolated from Minas Gerais State. 

In a brucellosis control and eradication program the 
use of an accurate surveillance and highly discrimin- 
atory typing method is essential to characterize an out- 
break and determine the source of infection and the 
transmission routes. The present results on typing mul- 
tiple B. abortus isolates from an outbreak originated 
from an outbreak in the same herd depicted the in vivo 
stability of the MLVA16 markers. The set of loci that 
comprise the MLVA16 demonstrated to be very stable, 
even when assessed over the time span of one year and 
four months, since the field strain mainly responsible by 
the outbreak (W) and the RB51 vaccine strain recovered 
from vaccinated animals showed unchanged MLVA16 
profiles. These findings are extremely important because 
they definitely confirm the ability of the MLVA16 to estab- 
lish correct epidemiological correlations, since, besides 
having a high discriminatory power, MLVA16 markers 
were also stable under natural selection pressure exerted 
by the host, during the sixteen months assessed. Thus, 
these results increase the confidence in the traceback 
established from the results of MLVA16 and further 
support this technique as the choice one for typing 
B. abortus. 

Futhermore, the examination of in vitro stability of the 
B. abortus RB51 vaccine strain, B. abortus strain 2308 
and B. abortus field isolates by serial passages in culture 
medium showed no change in MLVA16 profile [9,16]. 
Likewise, the analysis of B. abortus S19 vaccine strain 
from different batches of different manufacturers did not 
reveal significant differences in MLVA16 pattern [17]. 
Concerning the evaluation of in vivo genetic stability of 
MLVA16 loci, it has also been demonstrated that passage 
of B. abortus RB51 in cattle and of B. abortus 2308 in 
mouse did not lead to changes in any marker of 
MLVA16 [9]. Nevertheless, minor changes in VNTR pat- 
tern were observed in in vitro passages of B. abortus 544 
and in genotyping of multiple B. abortus isolates from the 
same outbreak [9]. Her et al. [9] found different allelic pro- 
files in seven of twenty-three herds in which more than 
one B. abortus isolates were obtained. Those different 
genotypes from same outbreak showed mutations only 
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in the loci Bruce 30 and 43, which did not seem to affect 
the identification of a possible common origin of the 
strains [9]. However, our data showing the high in vivo 
stability of the MLVA16 loci have as main findings over 
previous data the large numbers of isolates (37) ob- 
tained from the same source, the long period of time 
assessed (16 months) and the typing of strains coloniz- 
ing different sites (mammary gland or reproductive 
tract), and therefore under different selective environ- 
mental pressures. 

Conclusions 

The results of the present study on typing of multiple 
isolates from different clinical specimens originating 
from the same outbreak over a sixteen month period in- 
dicate the in vivo stability of the MLVA16 markers, a 
low genetic diversity among B. abortus isolates and the 
usefulness of MLVA16 for epidemiological studies of bo- 
vine brucellosis. 

Methods 

Outbreak description 

The experiment was conducted during an outbreak of 
brucellosis in a cattle herd located in Matozinhos, Minas 
Gerais, Brazil. The herd was composed of 705 Holstein 
dairy cows, reared in an intensive system. All animals in 
the herd were vaccinated with S19 between 3 to 8 month 
of age and tested for brucellosis twice a year according 
to PNCEBT [3] [Rose Bengal Plate Agglutination Test 
(RBPAT) as a screening test and the Standard Tube 
Agglutination Test (STAT) and 2-Mercaptoethanol Test 
(2ME) as a confirmatory tests]. During 2008, several lots 
of heifers were acquired in order to increase the milk 
production of the herd. Prior to introduction into the 
herd, all heifers were serologically tested [3] and only 
lots with brucellosis-negative animals were acquired. 

Suspicion of brucellosis in the herd started with an ab- 
normal high abortion rate (15%) observed by the end of 

2008. Brucellosis was confirmed serologically on March 

2009, and a control program based on mass vaccination 
with RB51 and a test-and-slaughter policy was initiated. 
PNCEBT demands the compulsory S19 vaccination of 
female calves aged 3 to 8 month and recommends the 
voluntary RB51 vaccination of heifers older than 8 month 
and cows in B. abortus infected herds. All female ani- 
mals over eight months of age were revaccinated with 
RB51, including pregnant heifers and cows. The herd 
was revaccinated with RB51 on two different occasions: 
the first in May 2009 and the second in December 2009. 
Moreover, all cattle older than 24 months were monthly 
tested for brucellosis (RBPAT as screening test and 
STAT and 2ME as a confirmatory test) and the positive 
animals were culled [3] until three consecutive negative 
herd results were obtained. 



The experimental protocol was approved by Ethics 
Committee on Animal Experimentation of Universidade 
Federal de Minas Gerais (CETEA/UFMG - Protocol 139/10). 

Clinical samples 

Three-hundred and seventy-five clinical samples were 
collected from 340 cows between January 2009 and May 
2010 during the brucellosis outbreak; those included 275 
vaginal swabs and 100 milk samples. Of thirty-five 
animals were collected both, vaginal swabs and milk 
samples. Vaginal swabs were collected in Stuart medium 
from all cows immediately after abortion (twenty-two 
samples) or delivery, whereas milk samples were obtained 
from all cows that had 60 days or less of the last partur- 
ition or abortion at 57 days after the first RB51 vaccin- 
ation (07/27/2009). For each animal, 50 mL of milk was 
collected after discarding of the first jet (milk from the 
four quarter was mixed). Both clinical samples were 
stored at -20°C until processing. 

Culture conditions 

Vaginal swabs were thawed and directly plated onto dupli- 
cate tryptose agar plates (Difco, USA) with antibiotics 
(Farrell's selective supplement) (Oxoid, UK) [31], and they 
were also inoculated in 10 mL of enrichment medium 
(tryptose broth supplemented with Farrell's selective sup- 
plement) [22]. Milk samples were thawed and centrifuged 
at 2500 x g, for 15 minutes. The intermediated phase was 
discarded and the supernatant was mixed with the pellet. 
Aliquots of 100 uL of each mixture was immediately inocu- 
lated onto duplicate tryptose agar plates (Difco, USA) with 
antibiotics (Farrell's selective supplement) (Oxoid, UK) [31]. 
Another aliquot of 1 mL of each mixture was diluted in 
9 mL of enrichment medium. The enrichment media inoc- 
ulated with vaginal swabs or milk were incubated at 37°C 
for seven days in 5% of C0 2 and then inoculated onto tryp- 
tose agar plates with antibiotics [22]. All plates were incu- 
bated in 5% of CO z at 37°C for at least 14 days [31]. 

Identification and biotyping of Brucella spp. isolates 

Brucella spp. solates were identified to genus based on 
colony morphology, positive tests for urease, catalase 
and nitrate reduction, and negative tests for motility and 
citrate utilization [31,32]. Classification into species and 
biovar was performed by H 2 S production, agglutination 
with sera anti-A, anti-M, anti-R and acriflavin, growth 
in C0 2 and 0 2 atmospheres and sensitivity to thionin 
(20 ug/mL and 40 ug/mL), basic fuchsin (20 ug/mL), 
and rifampicin (200 ug/mL) [31-33]. 

Suspected Brucella spp. colonies were also suspended in 
100 uL of Tris EDTA (10 mM Tris-HCl, 1 mM EDTA, 
pH 8.0), inactivated at 85°C for 2 hours and subjected to 
genomic DNA extraction [34]. DNA quality and concen- 
tration were determined by spectrophotometry [35]. 
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Polymerase chain reaction (PCR) for amplification of 
the gene bcsp31 [18] and AMOS -enhanced PCR [19] were 
also used for identification of isolates [see Additional 
file 1; Table SI]. Since B. abortus biovar 3 strains are 
not identified by AMOS -enhanced PCR, these strains 
were also analyzed by Bruce-ladder multiplex PCR [20] to 
confirm their identification as B. abortus [see Additional 
file 1: Table SI]. 

Reference strains, used as control in different procedures, 
were: B. abortus biovar 1 544 = ATCC 23448 T ; B. abortus 
biovar 2 ATCC 23449; B. abortus biovar 3 Tulya = ATCC 
23450; B. abortus biovar 4 292 = ATCC 23451; B. abortus 
biovar 5 3196 = ATCC 23452; B. abortus biovar 6 870 = 
ATCC 23453; B. abortus biovar 9 C68 = ATCC 23455; 
B. abortus biovar 1 S19; B. abortus biovar 1 RB51; B. 
melitensis biovar 1 16 M = ATCC 23456 T ; B. ovis Reo 
198; B. suis biovar 1 1330 = ATCC 23444; Escherichia coli = 
ATCC 25922; £ coli B41; Listeria monocytogenes = ATCC 
19115; Pseudomonas aeroginosa = ATCC 27853; Sal- 
monella enterica serovar Typhimurium = ATCC 14028; 
and Staphylococcus aureus = ATCC 29213. 

MLVA 16 genotyping 

DNA from each strain was genotyped by a subset of 16 
tandem repeat loci (MLVA16) as previously described [6,7] 
[see Additional file 1: Table SI]. The method analyses 16 
VNTR loci divided into three panels: panel 1 composed of 
eight minisatellites (Bruce06, Bruce08, Brucell, Brucel2, 
Bruce42, Bruce43, Bruce45 and Bruce55); panel 2A com- 
posed of three microsatellites (Brucel8, Brucel9 and 
Bruce21); and panel 2B with five microsatellites (Bruce04, 
Bruce07, Bruce09, Brucel6 and Bruce30) [6,7]. 

The amplified products were submitted to electrophor- 
esis in 2% or 3% agarose gel, for the mini and microsatel- 
lites, respectively, in Tris-borate-EDTA IX (TBE) buffer, 
stained with 0.5 mg/mL ethidium bromide, visualized 
under UV light, and photographed (ImageMaster VDS, 
Phamarcia Biotech, Sweden). DNA size markers 100 bp 
(100 bp DNA Ladder, New England Biolabs, USA) and 
25 bp (25 bp DNA Step Ladder, Promega, USA) were used 
to estimate the tandem repeat unit length. 

Each estimated band size was converted into number of 
repeat units for each locus, with aid of the software BioNu- 
merics 5.1 (Applied Maths, Belgium), as described [6]. Clus- 
tering analysis was performed using the same software based 
on the category coefficient and the unweighted pair group 
method with arithmetic mean (UPGMA) algorithm [6]. 

Additional file 



Additional file 1: Table SI. Primers used in the present study. 
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